Memory effects in adiabatic quantum pumping with parasitic nonlinear dynamics 
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The charge current adiabatically pumped through a mesoscopic region coupled to a classical 
variable obeying a nonlinear dynamics is studied within the scattering matrix approach. Due to 
the nonlinearity in the dynamics of the variable, an hysteretic behavior of the pumping current 
can be observed for specific characteristics of the pumping cycle. The steps needed to build a 
quantum pump working as a memory device are discussed together with a possible experimental 
implementation. 

PACS numbers: 73.23.-b,72.15.Qm 



I. INTRODUCTION 



The quantum pumping proposed by Thouless 1 is a 
phase coherence effect able to pump dc current by the 
out-of-phase adiabatic modulation of two system param- 
eters. The phase difference tp of the external signals pro- 
duces a charge current proportional to ui sin(</?) , uj be- 
ing the adiabatic pumping frequency. After the pioneer- 
ing work of Thouless, a scattering matrix approach to 
the adiabatic quantum pumping has been proposed by 
Brouwer—. Following the Brouwer formulation, the no- 
tion of non-interacting quantum pump 3,4 has been devel- 
oped in several directions. The most important applica- 
tion is the possibility to achieve a pure spin current injec- 
tor, which is an important goal in spintronics. Recently 
the idea of interacting quantum pump 5 beyond the free 
electrons model, has been introduced. To get read of non- 
trivial interaction effects, typically described in the con- 
text of the Anderson model^ for quantum-dot systems, 
a non-equilibrium Green's functions (NEGF) theory of 
the quantum pumping has been developed using several 
approximations. Alternatively to the NEGF approach, a 
powerful diagrammatic technique based on a generalized 
master equation method^ has been developed. More re- 
cently, another interesting study^ was performed aiming 
at generalizing Brouwer's formula for interacting systems 
to include inelastic scattering events. 
From the experimental 9 point of view pumping effects 
have been analyzed in confined nano-structures, as quan- 
tum dots, where the realization of a periodic time- 
dependent potential can be achieved by modulating gate 
voltages. 

Very recently it has been proposed the idea that the pres- 
ence of a classical state variable inside the scattering re- 
gion of the system can strongly affect the form of the 
pumped current. This effect is produced by the linear 
dynamics of a classical variable which modifies the phase 
relations introducing a dynamical phase shift 10 . 
In this work we study the charge current pumped through 
a mescoscopic region coupled to a parasitic classical vari- 
able governed by a weakly nonlinear dynamics supposed 
adiabatic on the typical time scales of the electron trans- 
port. 



The organization of the paper is the following: In SecllTl 
we summarize the scattering matrix approach to the 
quantum pumping as formulated by M. Moskalets and M. 
Biittiker in Ref.Q. In SeclIIIIwe derive a general expres- 
sion for the current pumped through a mesoscopic region 
containing a parasitic internal degrees of freedom gov- 
erned by weakly nonlinear dynamics. In Sec lIVI the spe- 
cific case of the Duffing nonlinear oscillator is discussed, 
while the conclusions are given in SeclVl 



II. SCATTERING THEORY OF QUANTUM 
PUMPING 



In the standard theory of the quantum pumping (i.e. 
without internal classical dynamics) a mesoscopic scat- 
terer coupled to two external leads (left/right) taken 
at the same temperature T and electrochemical poten- 
tial /x w Ep, Ep being the Fermi energy, is consid- 
ered. When the scatterer is subjected to two control- 
lable parameters Xi, whose adiabatic time dependence 
is Xi(t) = Xf + Xf sin(wt + ipi), the mesoscopic phase 
coherent sample is characterized by the scattering ma- 
trix S(Xx(t), Exploiting an instant scattering de- 
scription and assuming that the pumping amplitudes Xf 
are small quantities compared to Xf, the scattering ma- 
trix takes the following form: 



S(X 1 (t),X 2 (t))^S(Xf,X, 



(1) 



where s v are related to the parametric derivatives of the 
scattering matrix. Under equilibrium condition of the 
external leads, the electrons with energy E entering the 
scattering region (m-states) are described by the Fermi 
distribution p a n (E) = f (E) = (I + cxp(-^))-\ where 
a G {left, right}. On the other hand, the distribu- 
tion function f^ lt {E) of the outgoing particles leaving 
the mesoscopic region and entering the reservoir (namely 
crai-states) is affected by the interaction with the oscil- 
lating scatterer allowing the electrons to absorb or emit 
an energy quantum hu. The latter mechanism changes 
the initial distribution function fo{E) and produces a 
non-equilibrium distribution responsible for the genera- 
tion of a finite particles current. To compute the non- 
equilibrium distribution of the outgoing particles fa Ut (E) 
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one introduces two kinds of operators: The fermionic field 
operator a a which annihilates an incoming state in the 
lead a and the annihilation operator of the outgoing state 
b a in the same lead. The operators a a and b a written in 
energy representation are related to the scattering matrix 
by the relation 



Sr,,a/3ap(E + r)Hu), (2) 



where S° ap = S(X^,X°) afs . Exploiting Eq.© the non- 
equilibrium distribution f° ut (E) can be computed as 
(b\(E)b a (E)) , being (...) the quantum-statistical aver- 
age. According to this, we obtain 

f° ut (E)=Y,\ S ° a p\ 2 fo(E)+ KatfME + r,^), 

(3) 

where we have used the relation (a) a (E)di3(E')) = 
5 a /3S(E — E')fo(E). If we define as positive the current 
I a flowing from the scatterer to the lead a we can write 



1° = T 



dE[f^(E)-f (E)]. 



(4) 



The above formula can be explicitly evaluated in the 
small w limit and the expression of the current I a be- 
comes: 



In 
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2^^ 



{\s~,<xp\ 



(5) 



In case of a quantum pump working with two time mod- 
ulated parameters, the dc-current I a can be written in 
the following final form: 



In 
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( 6 ) 

When the quantum pumping is performed by TV pa- 
rameters (multi-parametric case) with time modulation 
X n (t) = X® + X% sin(o;i + (p n ), the charge current con- 
tains terms proportional to X^X^ sin(<y9 m — cp n ) showing 
the sensitiveness of the quantum pumping to all phase 
differences experienced by the carriers inside the scatter- 
ing region. 

The latter observation is very important in describing 
the case of an effective multi-parametric case. In this 
light when a two-parameters Thouless pumping is 
performed, the dynamics of a parasitic classical state 
variable describing some degree of freedom of the scat- 
terer can be excited by the external modulations. The 
effect of such parasitic state variable is the introduction 
of an internal (a priori unknown) dynamics forced 
by the pumping cycle and acting as a third pumping 
parameter (effective multi-parametric case). In analogy 
with the multi-parametric pumping discussed before, 
we thus expect additional contributions to the pumping 
current related to the parasitic dynamics. To be precise, 



in our language a parasitic degree of freedom is a 
classical state variable related to the internal dynamics 
(a priori unknown) of the scattering region. The notion 
of parasitic variable is conceptually similar to the one 
of parasitic impedance in electrical circuits theory: i.e. 
a degree of freedom not considered in designing the 
device whose dynamics may affect the performance of 
the system in an unexpected (usually unwanted) way. In 
that context, once the presence of an hidden impedance 
is recognized, an equivalent circuital model of the whole 
system can be obtained and a proper control of the 
device is achieved. In full analogy to this well known 
situation, we would like to build a scattering theory of 
the two-parameters Thouless pumping able to take into 
account the parasitic dynamics. The resulting theory 
is equally applicable when the internal dynamics is a 
priori known as is the case for a deformable quantum 
dot system as considered in Ref.[lO(. 
Thus, when a classical state variable is indirectly excited 
by the pumping parameters, an internal phase shift may 
activate additional contributions to the charge current. 
This mechanism becomes very interesting in the case of 
a weakly nonlinear dynamics able to produce non-trivial 
internal phase shifts. 



III. QUANTUM PUMPING WITH A 
PARASITIC NONLINEAR VARIABLE 

We now consider the case when the dynamics of the 
oscillating scatterer is affected by the internal classical 
dynamics of a parasitic variable y(t). The nature of the 
variable y, if not recognized and included in the pumping 
model, does not allow to precisely control the system even 
though the external pumping parameters Xi(t) act as a 
driving force on the parasitic dynamics. In this way the 
scattering matrix S(X\ (t) , X^ (t) ; y(t)) depends paramet- 
rically on the pumping parameters and on the parasitic 
variable whose dynamics is controlled by the pumping 
cycle in a way which is a priori unknown. When the 
equation of motion of y(t) is linear (or when it can be 
linearized around the working point), the dynamics of 
the parasitic variable in the presence of external param- 
eters can be described in terms of the classical Green's 
function x(t — t')> 



dt' X (t-t')F(X 1 (t'),X 2 (t% (7) 



F(X\, X-z) being the forcing term of the differential equa- 
tion. Once the dynamics has been formally solved the 
charge current can be computed as done in Eq.(j6|) or in 
the multi-parametric case. In the latter case, already dis- 
cussed in Ref.[l0|, the dynamical phase shift 4>d induced 
by the Fourier transform |x|e^ D of x(t ~ t') plays an 
important role in generating additional contributions to 
the pumping current. Such additional terms, weighted 
by \x\, can give information on the differential operator 
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related to the parasitic dynamics configuring the quan- 
tum pumping as a probe of internal dynamics. 
When the internal variable follows a nonlinear dynamics 
all the above considerations are still valid even though 
new features appear. We consider here the simplest non- 
linear dynamics in which the interaction between the 
electronic degrees of freedom and the parasitic variable 
is governed by the equation of motion of a forced Duffing 
oscillator— (in dimensionless units): 

y + Py + y + e{-ay + by :i + J2 9iX" sin(r+ &)] = 0, (8) 

i=l,2 

where the gt are gain functions, a is a detuning term able 
to modify the resonance frequency of the system, while 
b > is the nonlinear termi^. Furthermore, the time r 
is normalized to the resonance frequency fio of the har- 
monic oscillator (obtained for e = 0) while the pump- 
ing frequency uj is tuned on resonance (i.e. u) = VLq). 
In the weak interaction limit (e <C 1) one can measure 
the linear damping in terms of e making the substitu- 
tion ft — ¥ e/3. In this way a weakly nonlinear system 
is obtained which can be treated within the two-timing 
perturbation theory-^. Namely two time scales are intro- 
duced, one is the fast time scale t = r and the other is the 
slow time scale T = er, while the approximate solution 
of the problem is written as y(t, T) ps yo(t, T) + eyi(t, T). 
The derivative with respect to r, i.e. d/dr, has to 
be substituted by the operator dt + edr, while d 2 /dr 2 
(up to the first order in e) takes the form d 2 + 2ed t ,T- 
With the above rules the zero-th order of the Duffing 
equation reduces to dfyo + yo — 0, which is solved by 
y (t, T) = R(T) cos(t + $(T)). Using y Q (t, T) within the 
first order equation and imposing vanishing secular terms 
one obtains the following set of equations for R(T) and 
$(T): 



2R' + (3R = fi cos($) + f 2 sin($) 
3 

2&R + aR- -bR 3 = -fx sin($) + f 2 cos($), 



(9) 



where fx = g x X<( cos(</>) + g 2 X% and f 2 = giXf sm{tp) 
( here we fixed <p\ = and ip 2 — f), while R 1 and <&' 
represent the first derivatives with respect to T. The 
above equations completely characterize the zero-th or- 
der solution of the problem allowing to derive the pump- 
ing current produced when a weak nonlinear dynam- 
ics is present. When the pumping terms are absent 
(fi = H = 0) the stationary solution is y e = and 
the static scatterer condition is reached. When non- 
vanishing pumping terms are present the equilibrium 
condition R 1 = 0, $' = can be written as follows: 



R: 



a - ~;bR 2 e 
4 . 



fi 



(10) 



COs($ e 

sin($ e 



= M 



0Re 



-aR c + (3/4:)bRl 



where M = (I/x — icr y f 2 )/ y/ fx + /f is a rotation matrix, 
being a y a Pauli matrix. Once the equilibrium radius R e 
is known solving the first equation in (|10j) the dynami- 
cal phase shift $ e is determined by the second equation 
in (|10[) . In this way the asymptotic zero-th order solu- 
tion is completely determined, y e (j) = R e cos(r + $ e ). 
Using the asymptotic form of the scattering matrix 
S(Xx (t), X 2 (t); y e (r)) and proceeding as in the standard 
case, one obtains the formula of the pumped current: 



1-1° 

1 a — 1 a 



euj ■ 
2^ 



^{cos($ e )[ J R e X;A^cosfe)] + 

0,3 

+ sin($ e )[ J R e X;A^sin(^)]}, (11) 



where = ^{(dySapdxjS*^}, while 1° is the stan- 
dard pumping term as given in Eq. (|6]) . Eq. (jlll) represents 
the main result of this work. Its validity goes beyond the 
Duffing dynamics we are discussing here. Indeed every 
weakly nonlinear problem of the form 



x + x + eh(x, x, t) = 



(12) 



has a zero-th order asymptotic solution of the form 
x e (r) — R e cos(t + $ e ), being R e and $ e determined 
by the precise form of the nonlinear time-periodic func- 
tion h(x,x,r). Thus Eq. pTj) has to be considered valid 
for the whole class of weakly nonlinear problems. From 
the experimental point of view, the parasitic nonlin- 
ear dynamics of the variable y e {j) can be probed us- 
ing Eq. triTl) and considering the quantities R e sin($ e )A^ 

and R e cos($ e )A J a g as fitting parameters. For the fitting 
procedure it can be useful to estimate d y S a p by using 
the Fisher-Lee relation 



S a p {E F ,T) = 8 a0 -i\jT a (Xx ,X 2 )Tp {Xx , X 2 )G r (E F , r) , 

(13) 

where G t (Ef,t) is the instantaneous retarded Green's 
function associated to the scattering region assumed of 
resonant-like form 

G r (E F ,T) = [£; F - £ ( 2 ; e (r))+*^r s (X 1 (T),A 2 (T))/2]- 1 , 

while the tunneling amplitudes T s are related to the 
Green's function self-energy. In this way we obtain a 
model for the parametric derivative of the scattering ma- 
trix with respect to the parasitic variable: 



A,/ror°ex P (-2*e) 



(£f- £o ) 2 + (£ s iW 



(15) 



where A = {d y s(y)) , 6 = arctan(E. r°/2]/[£ F - e ]), 
£o being the energy level of the resonant state inside 
the scatterer in the absence of pumping terms. When 
Eq ps Ep and considering T° « T° R = T , one obtains 
the simple relation (d y S a p)o ~ i^/^a leading to the 
useful estimation A^ (A/r )3 f i{(i9x J '5a/3)o}- Within 
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FIG. 1: Equilibrium radius R e as a function of the detuning 
parameter a by setting the remaining parameters as follows: 
/3 = 0.75, b = 3.5, g±X? = g 2 X% = 0.35. The curves have 
been obtained from top to bottom by setting the pumping 
phase as ip — 0, 7r/4, n/2, 37r/4. 



the resonant hypothesis considered above the nonlinear 
dynamics affects the pumped current mainly via the 
terms i? e sin($ e ) and i? e cos($ e ). Thus in case of 
bistability of the forced nonlinear dynamics an hysteresis 
of the pumped current induced by R e is expected. In 
the latter case, the quantum pump displays a memory 
effect controlled by the features of the pumping cycle. 
To explain this point we now go back to the analysis 
of the forced Duffing oscillator governing in our specific 
example the dynamics of y. 



IV. RESULTS FOR THE DUFFING 
OSCILLATOR 




1 2 3 4 5 6 

<P 

FIG. 2: Critical value (full line) of the nonlinear term b c as 
a function of the pumping phase <p setting the remaining pa- 
rameters as follows: /3 = 0.75, b — 3.5, giXf = g 2 X 2 = 0.35. 
The dashed line represents the the actual value of the nonlin- 
ear term b = 3.5. When the full line lies below the dashed one 
the condition b > b c is fulfilled and an hysteretic behavior is 
detected in Fig[TJ For b < b c , R e (a) is a single valued func- 
tion. The vertical lines correspond to the pumping phases 
chosen in Fig[TJ 



that the phase space of the oscillator contains two sta- 
ble limit cycles characterized by different radius R^ and 

(2) 

Re , a third unstable limit cycle being present between 
the stables one. In this way the state of the parasitic 
variable can be controlled by using the detuning a (or 
equivalently the pumping frequency tu) and the pumping 
phase tp. Furthermore, the change of state of the parasitic 
variable produces an observable variation of the pumping 
current allowing us to modify and detect the state of the 
parasitic variable y. In this way a quantum-pump-based 
memory device is obtained. The critical value b c of the 
nonlinear term can be computed exactly differentiating 
the first relation in (ITU)) with respect to a and impos- 
ing diverging derivative for R e (a) (i.e. dR e /da — > ±oo). 
After some additional algebra, one obtains the following 
expression for b c (ip): 



The numerical analysis of first equation in (|10[) is 
shown in FigUJ where the equilibrium radius R e is re- 
ported as a function of the detuning parameter a and 
for different pumping phases <p. Starting from the up- 
per curve obtained for ip = we observe an hysteretic 
behavior as a function of a. Increasing the value of the 
pumping phase the hysteresis disappears when the value 
ip ps 7r/2 is reached. The hysteretic behavior of the Duff- 
ing oscillator is controlled by the nonlinear term b. In 
particular for values of b below the critical value b c the 
curves R e (a) are single valued functions, while when b c 
is exceeded (i.e. b > b c ) a bistable behavior is detected. 
Since in our case b = 3.5 is fixed, we have to assume that 
the pumping phase ip can vary b c . Very interestingly the 
above observation implies that the shape of the pump- 
ing cycle may affect directly the stability of the nonlinear 
parasitic variable. In particular the hysteresis indicates 



32V3/3 3 
27[A 2 + /I 



(16) 



where the dependence on the pumping phase ip is 
hidden in the forcing terms fi and f%. The critical 
curve b c (ip) is reported in Figf2] as a function of the 
pumping phase setting the remaining parameters as 
done in FigUJ The dashed line represents the actual 
value of the nonlinear term (i.e. b = 3.5). The analysis 
of the figure shows that when the critical curve b c (<p) lies 
below the dashed line (i.e. when b > b c ) an hysteretic 
behavior affects the system, while an increasing of the 
pumping phase beyond ip s» 7r/2 destroys the bistability. 
Thus we can manipulate the hysteretic threshold b c by 
modifying the features of the pumping cycle (i.e. Xf, ip). 
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FIG. 3: Hysteresis loop for the radius R e as a function of the 
detuning parameter a by setting the remaining parameters as 
follows: /3 = 0.75, b = 3.5, giX? = g 2 X£ = 0.35, tp = 0. 



A. A quantum-pump-based memory device 

In this Section we discuss how a memory device based 
on the quantum pumping works. In order to build a 
memory device a scattering region containing a nonlin- 
ear parasitic dynamics is needed. For instance, a para- 
sitic dynamics can be activated by an unknown capaci- 
tive coupling between the scatterer and a charge density 
located on the mesoscopic substrate which can be usu- 
ally modeled by a classical RLC circuit. Another more 
controllable way to produce nonlinear parasitic dynam- 
ics is provided by a scatterer able to elastically react in 
the presence of a charge density on it. Such system is 
a nanoresonator coupled to external leads whose experi- 
mental realization is reported in lief. If. while its use as 
deformable quantum dot in the quantum pumping con- 
text has been reported in Ref . [lfj. Following a derivation 
similar to the one given in Ref. [l0( and going beyond the 
linear terms in expanding the charge density, one ob- 
tains the center of mass dynamics of the scatterer (i.e. 
the parasitic degree of freedom dynamics) in the form 
given in Eq.©. Thus Eq.© has to be considered the 
minimal nonlinear model of the center of mass dynam- 
ics, while the term e[— ay + by 3 + J2i=i i 9iX£ sin(r + tpi)] 
has to be meant as the electrostatic force acting on the 
nanoresonator with the capacitive coupling provided by 
e. Once the center of mass dynamics has been recog- 
nized in the form of the Duffing oscillator all the argu- 
ments concerning its stability are the same as discussed 
in Sec lIVI In particular, we now study the case in which 
the system presents a bistable behavior which is needed 
to obtain a memory device. In Fig(3]we present the hys- 
teresis loop obtained setting the system parameters as 
follows: (3 = 0.75, b = 3.5, g x Xf = g 2 X$ = 0.35, tp = 0. 
It represents the equilibrium radius R e as a function of 
the detuning parameter a which can be experimentally 



controlled by using gate voltages or acting on the fre- 
quency of the pumping cycle. Starting from the point 
A of the loop and decreasing the value of a, we move 
along the lower branch until the point B is reached. As 
evident, from analyzing the second and third term in 
Eq. (fTTj) . along AB the system presents a continuous vari- 
ation of the pumping current since R e and $ e are con- 
tinuous function of a. Further decreasing a, moving from 
the point B, produces a jump toward the point C lo- 
cated on the upper branch (i.e. CD) of the hysteresis 
loop. This jump, corresponding to a sizeable change of 
i? e , can be detected by a strong variation in the pump- 
ing current if the second and third terms in Eg. (1111) arc 
of the same order of magnitude compared to the pump- 
ing term. Once the point C is reached, the loop can 
be closed by increasing a until the state D is achieved. 
Staring from D, a further increasing of the detuning a 
produces a big jump toward the initial state A of the 
hysteresis loop. The DA jump (as also the BC jump), 
detectable in the pumping current, indicates the change 
of state of the parasitic dynamics (e.g. the center of mass 
evolution of the nanoresonator) . In this way the bistable 
state of the parasitic variable is recorded by two refer- 
ence pumping currents taking different values pertaining 
to the upper or lower branch of the hysteresis loop. The 
above procedure constitutes the working principle of a 
memory device based on quantum pumping. 



V. CONCLUSIONS 

In conclusion, we have shown an unconventional ap- 
plication of the quantum pumping in which a parasitic 
nonlinear variable affects the pumping cycle introducing 
interesting memory effects controlled by the critical pa- 
rameter b c . The hysteretic threshold b c can be manipu- 
lated acting on the shape of the pumping cycle and the 
state of the parasitic variable y can be detected measur- 
ing the pumping current. The features above are useful in 
obtaining a memory device based on quantum pumping 
which could be realized using a nano-resonator similar 
to the one reported in Ref. 14]. Apart from the tech- 
nological application of our study, nanoelectromechani- 
cal systems are ideal candidates to validate the present 
theory and explore the parasitic dynamics. Studies of 
nanoelectromechanical systems are usually performed by 
reducing their dynamics from that of a continuous model 
to that of the center of mass dynamics, the so called "res- 
onator" . However in the real situation the resonator is an 
extended object whose center of mass motion is affected 
by different vibrational modes of the structure. Thus, 
due to the complexity of the problem is not possible to 
fix a priori the parameters of the reduced model (e.g. 
the damping coefficient, the resonance frequency, etc) or 
recognize the structure of the nonlinear contributions to 
the damping term. In fact the total energy is differently 
stored among the vibrational degrees of freedom depend- 
ing on the working point of the resonator. Thus, even 
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though the resonator dynamics can be taken into account 
in designing the device, to fix the time evolution of the 
mechanical degrees of freedom a precise knowledge of the 
capacitive coupling as a function of the resonator posi- 
tion is needed. Our proposed method can be useful in 
identifying the model for the capacitive coupling and of- 
fer a method to probe the parasitic nonlinear dynamics. 
A further example of parasitic dynamics is constituted by 
the capacitive coupling of a classical background charge 
to the scattering region. Such situation can be easily 
understood as the capacitive coupling between the scat- 



tering region and an effective RLC circuit. However, also 
in this simple case, the parasitic dynamics can be imag- 
ined as the linearization of some nonlinear model around 
the working point whose dynamical parameters (capaci- 
tance, differential resistance, etc) can be tested by using 
our theory. 

Finally, a variety of other useful applications such as 
frequency synchronization, frequency mixing and con- 
version, and parametric amplification, could be imple- 
mented by our proposal playing an important role for 
nanomechanical devices. 
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